Search for Excited Quarks in qq — > 77 at the LHC 
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If quarks are composite particles, then excited states are expected to play a role in the Large 
Hadron Collider phenomena. Concentrating on virtual effects, and using a large part of the CMS 
detection criteria, we present here a realistic examination of their effect in diphoton production at 
the LHC. For various luminosities, we present the 99% confidence limit (CL) achievable in A — M qt 
parameter space where A is the compositeness scale and M q * the mass of the state. For a q* of 
mass 0.5 TeV, A < 1.55 (2.95) can be excluded at 99% CL with 30 (200)fb _1 integrated luminosity. 

PACS numbers: 12.60.Rc, 13.40.-f, 13.85.Qk 



I. INTRODUCTION 

The replication of fermion families alongwith the mass 
hierarchies and mixings has led one to speculate about 
the possibility of quark-lepton compositeness, namely 
that the Standard Model (SM) fermions are not elemen- 
tary at all. The fundamental matter constituents in such 
theories, very often termed preo7is[l|, experience an hith- 
erto unknown force on account of an asymptotically free 
but confining gauge interaction 0, which would become 
very strong at a characteristic scale A, thereby leading 
to bound states (composites) which are to be identified 
as quarks and leptons. In most such models 0,13], quarks 
and leptons share at least some common constituents. 

If this hypothesis were to be true, it is possible, in- 
deed probable, that excited states of fermions exist at 
a mass scale comparable to the dynamics of the new 
binding force. In the simplest phenomenological mod- 
els 0, the excited fermions are assumed to have both 
spin and isospin 1/2 and to have both their left- and 
right-handed components in weak isodoublets (i.e. they 
are vector- like). Since these interact with the SM parti- 
cles, they may be produced at high-energy colliders and 
would decay back, radiatively, into an ordinary fermion 
and a gauge boson (photon, W, Z or gluon). Pair pro- 
duction of charged excited fermions could proceed via 
s-channel (7 and Z) exchanges in e + e~ collisions, while 
for excited neutrinos only Z exchange contributes. Al- 
though t-channel diagrams are also possible, they gen- 
erally give a negligible contribution to the overall pair 
production cross-section on account of the smallness of 
the cross-couplings [5]. However, this very same interac- 
tion between the excited state, its SM counterpart and a 
gauge boson may be used to singly produce such states 
(through both s- and f-channel diagrams). The four LEP 
collaborations have used these (and other) modes to es- 
sentially rule out such excitations almost upto the kine- 
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matically allowed range [|. At the hera, on the other 
hand, both excited leptons and quarks may be produced 
singly through i-channel diagrams and these processes 
have been looked at without any positive results 

At the Tevatron, one may either pair-produce the ex- 
cited quarks (primarily through gauge couplings) or pro- 
duce them singly via quark-gluon fusion, provided the 
q*qg coupling strength is significant. A striking signal 
of the latter would be an enhancement in the dijet pro- 
duction rate with a peak in the invariant-mass distribu- 
tion. Whereas the D0 collaboration has also excluded 
the mass region 200 GeV < M q * < 720 GeV for excited 
quarks decaying to two jets Q, the CDF collaboration 
considered a multitude of decay channels, thereby ex- 
cluding the mass range of 80 GeV < M g * < 570 GeV 

The presence of such particles would change the phe- 
nomenology even if they were too heavy to be produced. 
Since the confining force mediates interactions between 
the constituents, it stands to reason that these, in turn, 
would lead to interactions between quarks and leptons 
that go beyond those existing within the SM. Well be- 
low the scale A, such interactions would likely be mani- 
fested through an effective four fermion contact interac- 
tion [ll], [l2j term that is an invariant under the SM gauge 
group. The D0 and the CDF experiments at the Teva- 
tron have searched extensively for excited quarks decay- 
ing to different final states as predicted by various mod- 
els, with the negative results translating to lower bounds 
on compositeness scale A. The D0 collaboration has put 
a lower bound of A > 2.0 TeV at 95% CL from an anal- 
ysis of dijet production The CDF collaboration has 
also put a lower limit of A > 2.81 GeV at 95% CL study- 
ing the qq — > ev process [HJ]. From a phenomenological 
study of flavor independent contact interaction for the 
diphoton final state, the lower bound for the LHC has 
been estimated to be A± > 2.88 (3.24) TeV at 95% CL 
for an integrated luminosity of 100 (200) fb' 1 

As can be readily appreciated, the different production 
modes (and decay channels, wherever applicable) probe 
different aspects of the effective theory that governs the 
low energy interactions of these excited states. In this pa- 
per, we seek to concentrate on one such property namely 
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the trilinear coupling of the excited quark to its SM 
counterpart and the photon. To be more precise, rather 
than seeking to actually produce these excited states, we 
would like to investigate their role in photon pair produc- 
tion at the LHC. Analogous to the process e+e - — > 77(7) 
used to probe compositeness at LEP, such an exercise 
would complement the excited quark direct searches for 
the mass region above the kincmatical threshold. Since 
diphoton production is both a very simple final state and 
likely to be well-studied at the LHC, it is of interest to 
see how well can this mode probe compositeness. 

The rest of the paper is organized as follows. In the 
next section, we discuss the effective Lagrangian for the 
theory under consideration and the new physics contri- 
bution to diphoton production. In section III we dis- 
cuss various SM backgrounds for the signal. In sections 
IV and V respectively, we describe the event generation 
and photon candidate reconstruction. Isolation study for 
photon is discussed in section VI. Confidence limit calcu- 
lations and results are presented in sections VII and VIII 
respectively. The systematics is discussed in section IX, 
and in the last section we summarize this analysis with 
our conclusions. 



is not unexpected in an effective theory as the term in 
cq. ([1]) is only the first term and the loss of unitarity, to a 
given order, is presumably cured once suitable higher di- 
mensional operators are included. An equivalent way to 
achieve the same goal is to consider the Cj to be form fac- 
tors rather than constants. To this end, we shall consider 
the q*q"f vertex to be given by 



3*37*i(p) 



A V A 2 



V P 



(2) 



where Q denotes a relevant momentum transfer. It can 
be checked that, for Q 2 = s, unitarity is restored as long 
as the constant n > 1. In the rest of our analysis, we 
shall confine ourselves to a discussion of n = 1. While 
this might seem to be an optimistic choice, it is not quite 
so. As can be readily appreciated, such a form factor 
plays a non- negligible role only when Q 2 ~ A 2 . Since, 
at the LHC, we shall prove to be sensitive to A of the 
order of a few TeVs, clearly the form factor plays only 
a marginal role in the determination of the sensitivity 
reach. 



II. EXCITED QUARK CONTRIBUTION TO 
DIPHOTON PRODUCTION 



As our interest is not in the production of the excited 
states, but rather on their contribution to the diphoton 
rates at a hadronic collider, it suffices to consider only 
the relevant parts of the Lagrangian, namely the mag- 
netic transition between ordinary and excited states. In 
general, it is often parametrized by 
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where the index i runs over the three SM gauge groups, 
viz. SU(3), SU(2) and U(l) and g u Gf^ and Tf are 
the corresponding gauge couplings, field strength ten- 
sors and generators respectively. The dimensionless con- 
stants Ci are, a priori, unknown and presumably of order 
unity. Clearly, the phenomenology would depend con- 
siderably on the ratios of the constants Cj. For exam- 
ple, electromagnetic couplings (and hence such decays) 
of such fermions are forbidden if ci — e/ci. Thus, the 
search strategies would depend crucially on the strengths 
of these couplings. 

A further point needs to be noted here. In the event of 
any one of the c' { s dominating the others, the cross sec- 
tion for any process governed by the Lagrangian above 
would scale as some power of the ratio Cj/A. Thus, in 
such a case, it makes sense to eliminate Ci altogether in 
favour of the unknown scale A. Furthermore, with the 
Lagrangian of eq.fl]) being a higher dimensional operator, 
the cross sections would typically grow with the center 
of mass energy, consequently violating unitarity. This 
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FIG. 1: Production process for diphoton final state (a) Ex- 
cited quark mediated (b) SM production. 

With the introduction of the new (flavour-diagonal) 
vertex as in eq.©, the process qq — > 77 acquires a new 
contribution as portrayed in FigfT] The differential cross 
section for the partonic process now reads 
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where the SM result is recovered in the limit A — > 00. 
The new physics contribution to the differential cross sec- 
tion thus depends on only two parameters, namely A and 
the mass of the excited state M q * . For simplicity, we as- 
sume these to be flavour-independent (within a genera- 
tion, it obviously has to be so). For eq.fl]) to make sense 
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as an effective Lagrangian, the masses have to be less 
than A (Ref.[H[ requires that M r < A/Vfy- 
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FIG. 2: The contribution of new physics to the total cross 
section. 

In Fig. [21 we present the additional contribution to 
the total diphoton cross section accruing from the new 
physics terms in eq.(3). Note that, unlike the QED 
contribution, this additional contribution does not suf- 
fer from collinear singularities. 

Contrary to the case of actual production (and subse- 
quent decay) of the excited state, the case under consid- 
eration is not associated with any resonant peak. Nev- 
ertheless, the presence of the new contribution preferen- 
tially enriches the large invariant mass end of the dipho- 
ton spectrum. The exchange of a large mass particle in 
the t— and it-channels naturally enhances the high-pr 
photon sample. To improve the signal to noise ratio, we 
must then concentrate on such a phase-space restricted 
subset of the final state. 

As can be gleaned from a cursory examination of 
eq.(3), the aforementioned dependence of the new con- 
tribution on the photon px is not as extreme as that 
for the QED contribution. Thus, the imposition of such 
cuts as we will discuss later would not drastically change 
the shape of the iso-cross section contours as depicted 
in Fig [2j Consequently, the exclusion contours that we 
shall finally obtain would bear considerable similarity 
with those in Fig [21 



III. BACKGROUND 

Standard Model processes, understandably, produce a 
large background to the diphoton final state. The back- 
ground can be divided into two categories: 



• where two prompt photons are produced in the 
(hard) subprocess itself, and 

• in a 7 + jet sample, a jet with a large electromag- 
netic fraction (e.g, ir a ,uj,rj etc.) fakes a photon or 
a hard photon is produced in the process of frag- 
mentation. 

The first category is dominated by the Born-level pro- 
cess qq — > 77. An additional source of the diphoton final 
state is provided by the gg — ► 77 process induced by a 
box diagram. Although the cross-section for this pro- 
cess is relatively small compared to the Born production 
(in fact, much smaller if very forward photons were to 
be included) the much larger gg luminosity at the LHC 
energies implies that gg — > 77 can be quite important. 
Indeed, even after imposing our selection criteria (to be 
discussed later) of moderately low rapidities and high 
transverse momenta for the photons, the gg-initiated con- 
tribution is approximately 6.8% of the Born contribution 
(see Table I). 



TABLE I: Various SM cross-sections for P T >190 GeV and 
|r/| <2.7 at -Js = 14 TeV. P T , the CKIN(3) parameter in 
PYTHIA, is the Pt of the outgoing partons in center of mo- 
mentum frame in a 2 — > 2 hard scattering process. 



Process 


Cross-Section (fb) 


7 + 3 et 


48970 


qq — ► 77 (Born) 


76.05 


gg -> 77 (Box) 


5.18 



Apart from the Born and box processes, single photon 
production processes qg — > 7(7, qq — ► jg and gg — ► 75 
where a jet fakes a photon can be a major source of back- 
ground. We have considered all these processes for the 
background estimation. Although the probabilty of a jet 
faking a photon is ^ 10~ 3 — 10~ 4 , the cross section for the 
first two of these hard processes (qg ^79 ,qq — ► 73) are 
larger by a typical factor of 0(a s /a) apart from a sin- 
gle ratio of gluon to quark densities, thereby partly rec- 
ompensing for this suppression. The third process, viz. 
gg — ► 73, is once again box- mediated and significantly 
smaller than the other two. Similar considerations hold 
for the background from dijet production with both jets 
being identified as photons. While the dijet cross section 
is very large, isolation requirements reduce it drastically. 
Even a simple estimate, without a full simulation, shows 
it to be quite unimportant for the physics under investi- 
gation. 



IV. MONTE CARLO SIMULATION & CUTS 

To generate the signal as well as the background 
events, we have used the PYTHIA [TtJ event generator 
wherein the signal matrix element of Eq.([3]) had been 
properly incorporated inside the pythia framework. It 
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was also counterchecked with a parton-level Monte Carlo 
generator. We have used the CTEQ5L parton distribu- 
tions [HI, with a choice of Q 2 — s for the factorization 
scale. While generating events, the multi parton interac- 
tion (MPI), initial state radiation (ISR) and final state 
radiation (FSR) switches in pythia were kept "ON" . 




A (TeV) 



FIG. 3: Variation of cross section with A at ^/s — 14 TeV 

In view of the fact that the signal events preferen- 
tially populate the large transverse momentum part of 
the phase space, events were generated with Pr > 190 
GeV (CKIN(3) parameter) and |?7| < 2.7 respectively. 
This also rids us of a very large fraction of the SM events 
which, understandably, are peaked at small angles to 
the beam-pipe as also small transverse momenta. Fig. [3] 
shows the variation in cross section with A for a fixed 
value of M qif =1.0 TeV. Clearly, the variation is well- 
approximated by a A -2 contribution superimposed upon 
a constant (the SM value). This is reflective of the fact 
that, for large A, the new physics contribution is domi- 
nated by the interference term in Eq.([3]) rather than the 
pure A -4 term . Only if we had imposed harder cuts 
on the photons, would the latter term have dominated 
(albeit at the cost of reducing event numbers and hence 
the sensitivity). 

It must be noted at this stage that, in the final selec- 
tion, we have used the fiducial volume of the electromag- 
netic calorimeter of the CMS detector i.e. |?7| < 2.5 with 
1.444 < |?7| < 1.566 excluded on account of the insensitive 
region between the barrel and the endcaps[19|. 



V. PHOTON CANDIDATE 

Since the SM 7 + jet and jet-jet production processes 
form a significant background to qq — > 77 via q* ex- 
change, it is very important to understand the mecha- 
nism of a jet faking a photon. The identification of a re- 
constructed object as a photon candidate depends on the 
specific design of the detector and the reconstruction al- 
gorithm. Taking this into consideration, at the generator 
level, we have used a clustering algorithm to account for 



fake photons arising from jets [2(j. The CMS experiment 
uses PbWOi crystals for the electromagnetic calorimeter 
(ECAL). Each crystal measures about 22 x 22 mm 2 [2l| 
and covers 0.0175 x 0.0175 (1°) in the Arj - A0 space 
((/> being the azimuthal angle). For photon reconstruc- 
tion, we have used the "hybrid" algorithm [T§] . The first 
step is to find a seed above a certain minimum tranverse 
momentum threshold P^ ln of 5 GeV[2l|. Only electro- 
magnetic objects, i.e., 7,e + and e~ are chosen as seed. 
Subsequently, one looks for all electromagnetic particles 
around the seed in the r/ — <j) space where Arj and A<j) dis- 
tance from the seed object is at most 0.09. This extension 
is equivalent to 10 x 10 crystal size in the CMS detector. 
The CMS experiment uses 5x5 crystal size to form an 
energy cluster and nearby non-overlapping clusters are 
merged to reconstruct a photon candidate. However, in 
our effort to mimic this reconstruction process at the gen- 
erator level, we choose to be conservative and use only 
a 10 x 10 crystal. We define the momentum of a pho- 
ton candidate to be the vector sum of the momenta of 
the electromagnetic objects in such a crystal. A photon 
candidate will be either a direct photon or other electro- 
magnetic obejets such as tt° — ► 77, p° — > 77 etc. Events 
where the two highest E t photons have cos(# 7 i 7 2) > 0.9 
with # 7 i 7 2 being the opening angle between the two pho- 
tons, are not considered because they could merge into 
a single energy cluster in the real detector. We have 
compared our results with the fast detector simulation 
(FAMOS [11]) used for CMS experiment and they are 
found to be in good agreement. With this algorithm 
and requiring the photon to be isolated (to be discussed 
later), the estimated probability of a jet faking a photon 
in 7+ jet channel is ~ 10~ 3 — 10 -4 . The major sources of 
fake photons are tt° (~ 81%), rj (~ 12%) and lo (~ 3%), 
with only a small fraction coming from other sources. 



VI. ISOLATION VARIABLES 



In a detector, a photon is recognised as a local depo- 
sition of electromagnetic energy in a limited region in 
the rj-<p phase space. In practice, it is defined as elec- 
tromagnetic energy contained in a cone of a given size 
R = y A(j) 2 + Ar/ 2 with no associated tracks. Fake pho- 
ton signals arising from a jet can be rejected by requiring 
either the absence of charged tracks above a certain min- 
imum transverse momentum (P^min) associated with the 
photon or the absence of additional energetic particles 
in an annular cone (Ri SO ) around the photon candidate. 
We have considered two variables for the isolation pur- 
pose (a) the number of tracks (N tr k) inside a cone around 
the photon and (b) the scalar sum of transverse energy 
(Etsum) inside a cone around the photon. 
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A. Track Isolation 

We have considered "stable" charged particles e.g. 
7T , K ± , e ± and P ± as tracks. Of these, ir^ alone con- 
tribute ~ 80% of the total charged tracks. The con- 
tributions from stable charged particles other than the 
ones mentioned above are negligible. The distributions of 
the number of charged tracks with a requirement on the 
transverse momentum of the tracks pointing to either the 
leading photon or the second leading photon candidate 
and within a corresponding cone of size 0.35 are shown 
in Fig.HJ In the signal sample (although we demonstrate 
for a particular value of the parameters, the features are 
generic), both photon candidates are true photons and 
hence the distribution falls very rapidly. The situation is 
markedly different for the background. For a true 7+jet 
event, the second leading photon is usually the fake one 
and has a large amount of hadronic activity around it. 
Consequently, the distribution (in Fig. reaches a max- 
imum around 5-6 tracks and then falls slowly. To under- 
stand the shape of the background distribution in Fig.|4jz, 
it should be realized that a small fraction of such events 
would actually have the fake photon as the leading one. 
Since such photons have a large number of tracks around 
them, an extended tail as seen in Fig. \Qi results. The 
same effect leads to the rise in the background distribu- 
tion for the second- leading photon for Ntrk < 1 (Fig.©). 



In pp collisions at the LHC, one expects to have a large 
number of soft tracks from associated minimum bias and 
underlying events. The major sources of tracks in the 
case of a true photon case are ISR, FSR and MPI, while 
the \ow-P T rk (< 1.5 GeV) tracks emanate mainly from 
the debris of the colliding protons. If these tracks are 
counted, a true isolated photon emitted from a hard pp 
collision may also appear non-isolated, thereby reduc- 
ing the signal efficiency. To avoid such possibilities, soft 
tracks are cleaned up by requiring the tracks to have a 
Pt above a certain minimum threshold (-PjC?£n) - ^ n var i~ 
ous CMS studies P^mm typically varies between 1-2 GeV 

sum- 
in this analysis, we have considered several choices for 
P T % n , namely 0.0, 1.0, 2.0 and 3.0 GeV respectively, 
and for different isolation cone sizes. The signal efficiency 
and the signal over background (S/B) ratio were calcu- 
lated with these choices for P T r J? lin and for various N tr k 
possibilities. The results, for the second leading photon, 
are displayed in Fig. [5l As one can observe, for N tr k — 
0, as P$£ in is increased from 1.0 GeV to 3.0 GeV, the 
signal efficiency increases by more than 15% with only 
a small reduction in the S/B ratio. Although, allowing 
more tracks in a given cone size leads to an increase in 
the signal efficiency, the S/B ratio decreases drastically 
(see Fig. 15). 
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FIG. 4: Number of tracks for the signal and the background 
events with P^ rk > 3.0 GeV pointing (a) leading photon and 
(b) second leading photon candidates in a cone of size 0.35. 



FIG. 5: Effect of the minimum threshold for track Pt on the 
S/B vs efficiency distribution for the second leading photon. 

Understandably, neither the SM diphoton contribution 
(whether the Born or the box-mediated processes) nor 
the new physics contribution to the same are affected by 
the requirement of Nt r k—0- Only the 7 + jet background 
suffers. Fig(S] shows the corresponding distribution in 
Pt for the highest transverse momentum track emanat- 
ing from the second leading photon. Both the distribu- 
tions (signal and background) have been normalized to 
unity. Clearly, the background dominates the signal for 
^Tmin > 3.5 GeV, thus pointing out a means to reject 
a large fraction of the 7 + jet background. Only those 
events are accepted where neither of the photons have an 
associated track with Pt >3.0 GeV within the respec- 
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tive isolation cones (i.e. N tr k—0 for P T rk >3.0 GeV). 
Only the highest Pt track is considered because consid- 
ering lower Pt tracks may affect signal efficiency. Since 
this study has been done at the generator level we have 
chosen P T r £ in >3.0 GeV. 
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FIG. 6: Highest track Pt around the second leading photon 
for both signal and 7 + jet background. An isolation cone of 
size 0.35 has been used. 
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FIG. 7: Etsum for the signal and the background events 
around (a)the leading and (b)the next leading photons. 



Defined as the cluster of energy inside a cone Ai? from 
which the energy of the photon is subtracted, the vari- 
able Etsum can be used to discriminate against an event 
wherein a jet fakes a photon. Although, in a real detector, 
Etsum is separately accounted for in the electromag- 
netic and the hadronic calorimeters, due to limitations 
of a generator level study, we use a combined Etsum 
which is the scalar sum of transverse energy of the elec- 
tromagnetic and hadronic particles around the photon 
candidate. 

Fig. [7] shows the normalized Etsum distributions for 
the signal and the backgrounds. The main aim of this 
study is to optimize the Etsum isolation variable so as to 
reduce the background from 7 + jet events. The leading 
photons, expectedly, have similar distribution for the sig- 
nal and the background. For the second photon though, 
the behaviours are very different. Most of the 7 + jet 
events have Etsum > 5 GeV and by Etsum ~ 10 GeV, 
the S/B ratio is miniscule. 

In Fig. [51 we show the variation in signal efficiency 
and the S /B ratio for different cone sizes around the sec- 
ond leading photon. Each point corresponds to a dif- 
ferent Etsum threshold, varied in steps of 1 GeV be- 
ginning with 1.0 GeV. The final choice of the cone size 
and the Etsum threshold depends on the track isolation 
efficiency, the signal efficiency, and the S/B ratio. 




1- 

0.44 0.445 0.45 0.455 0.46 0.465 0.47 0.475 
S/B 

FIG. 8: Signal efficiency vs. S/B ratio for different cone sizes 
for different choices of the Etsum threshold around the sec- 
ond leading photon. 



C. Combined Isolation 

In Table II, we show various combinations of isolation 
variables for two different cone sizes. Since we aim to 
observe an excess of diphoton production over the SM 
expectations, it is rather important to have a large sig- 
nal efficiency. We have performed this study for a large 
number of A — M q * points for which the cross section is 
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TABLE II: Fraction of events surviving for signal and back- 
ground after applying isolation cuts on both photons (and the 
& rp criteria). Also shown is the S/B ratio. 



XL) en Tvtrb 


rpmax 
1 oU M 


rytrk 

Tmin 

(GeV) 


ga 
(%) 


Born 

(%) 


Box 

(%) 


7 + Jet 
(%) 


S a /B 


0.35 


4.0 


1.5 


75.53 


75.45 


71.86 


0.81 


0.828 






2.0 


80.52 


80.40 


76.84 


0.90 


0.824 






3.0 


83.33 


83.19 


79.57 


0.96 


0.821 




5.0 


1.5 


77.10 


77.05 


73.60 


0.86 


0.824 






2.0 


83.15 


83.05 


79.68 


0.98 


0.818 






3.0 


87.18 


87.19 


83.79 


1.09 


0.810 


0.30 


4.0 


1.5 


81.20 


80.99 


77.97 


0.97 


0.817 






2.0 


85.73 


85.59 


82.55 


1.07 


0.811 






3.0 


88.49 


88.32 


85.44 


1.15 


0.806 




5.0 


1.5 


82.25 


82.17 


79.10 


1.01 


0.813 






2.0 


87.48 


87.45 


85.45 


1.14 


0.805 






3.0 


91.30 


91.24 


88.39 


1.26 


0.798 



Here A = M qt = 2.5 TeV 



TABLE III: The number of events surviving for signal and 
background for L ln t = lfb -1 after applying the final selection 
criteria. 



R,so Ntrk E¥sum Prmin S a Born Box 7 + Jet Tot.Back. 
(GeV) (GeV) 



0.35 


4.0 


1.5 


46.32 35.23 


1.90 


5.67 


42.81 






2.0 


49.31 37.54 


2.03 


6.29 


45.87 






3.0 


51.09 38.85 


2.10 


6.70 


47.66 




5.0 


1.5 


47.33 35.98 


1.94 


5.99 


43.92 






2.0 


50.94 38.78 


2.11 


6.84 


47.73 






3.0 


53.54 40.71 


2.21 


7.56 


50.49 


0.30 


4.0 


1.5 


49.83 37.82 


2.06 


6.74 


46.62 






2.0 


52.55 39.96 


2.18 


7.46 


49.62 






3.0 


54.35 41.24 


2.26 


8.00 


51.51 




5.0 


1.5 


50.57 38.35 


2.09 


7.01 


47.46 






2.0 


53.67 40.83 


2.23 


7.93 


51.00 






3.0 


56.10 42.60 


2.34 


8.78 


53.75 


"Here A = 


1.0 TeV and M q 


» = 0.5 TeV 









slightly larger than qq — > 77 production cross section, or 
in other words those points for which there will be only 
a small excess over the SM background. Although we 
have used a simple approach, it is possible to have other 
criteria to select analysis points for the choice of final 
selection cuts. Based on the studies detailed above, the 
final selection cuts are as follows: 

• P} 1 > 200 GeV, P'rf 2 > 200 GeV; 



• l^'T 2 ! < 2.5 

• COS(P-yl-y2 ) < 0.9 
jtrk 



7 i, 7 2| 



g [1.4442,1.5666]; 



• N trk = for P|T fe > 3.0 GeV within R lso < 0.35; 

• Etsum < 5.0 GeV within R lso < 0.35. 

After the application of the fiducial volume and photon 
Pt criteria, the requirement on angular separation be- 
tween the photons removes only ~ 1% events. 

Table III shows the number of events surviving for sig- 
nal, Born, box, 7 + jet and total background for 1/6 -1 
of integrated luminosity after applying the final selection 
criteria. 

Fig. [9]shows the distributions for some of the variables 
for the generated signal and background events after the 
selection requirements are imposed. In Fig[Sf, 9* is the 
angle between the direction of boost of the diphoton sys- 
tem and each photon in the diphoton rest frame. 



VII. CONFIDENCE LEVEL CALCULATION 

As the q* appears only in the ^-channel, no resonance 
peak appears in the diphoton invariant mass distribu- 
tion. Rather, a discovery needs to be made from an ob- 
servation of enhanced rate in the diphoton channel as 
well as differences in the shape of diverse phase space 



distributions. In this analysis, we primarily use the in- 
formation contained in the invariant mass distribution 
to distinguish between two hypotheses, namely the sig- 
nal + background hypothesis (S+B) and the background 
only (B) hypothesis. We adopt a frequentist approach 
to determine the confidence level of rejecting the S+B 
hypothesis (the exclusion CL) in the absence of a signal. 
The histograms shown in FigJ^e) are used to generate 
two sets of Gedankenexperiments. To do so, we assume 
that the content of each bin in the histograms is Pois- 
son distributed. For every bin, we generate a Poisson 
random number, taking the original bin content as the 
mean. These Poisson fluctuated random numbers now 
represent the bin contents of a histogram which we call 
the outcome of a single Gedankenexperimcnt. One mil- 
lion such Gedankenexperiments are generated from the 
S+B histogram (and, similarly, from the B histogram). 
From each of these experiments we calculated our test 
statistic, namely 



Xs+b 



E 

j=i 



(dj - (S + B) t f 



(4) 



(and similarly for Xb)- Here, di is the number of events 
in the i th bin of the M 77 distribution as generated in a 
particular Gedankenexperiment and (S + B)i is the num- 
ber of events in the original histogram of Af 77 obtained 
from PYTHIA. The distribution of x 2 shows how the 
test statistic will be distributed over many repeated ob- 
servations of the mass histogram. In Fig. [TU1 the solid 
histogram shows the expected distribution of \ 2 if the- 
S+B hypothesis is true while the dotted one shows the 
X 2 distribution if the S+B hypothesis is not true. The 
most probable value of x 2 if S+B is false is given by the 
peak of the Xb distribution. The area, a of the x's+b 
curve to the right of this value is the probability of see- 
ing a x 2 value > x% (peak) if the S+B hypothesis is 
true. For every point in the (A, M q *) plane satisfying 
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1 - a > 99%, the point is rejected at 99% CL. 

In calculating the x 2 j on ly bins with large significance 
are used. These have large bin contents and the latter 
can be safely assumed to be Gaussian distributed. As a 
consequence, the x 2 statistic detailed above is equivalent 
to a log likelihood statistic for this analysis. 

Since we have used histograms generated from 
PYTHIA as our input for the CL extraction there is sta- 
tistical uncertainty associated with the procedure, i.e., in 
a repeat of this MC study the position of the x\ P e& k 
will fluctuate, resulting in a different value of a. However 
at 1 — a = 99%, this fluctuation is estimated to be less 
than 0.5% on either side of the peak. 



VIII. RESULTS 

Fig. [TT1 shows the A — M q * parameter space which can 
be excluded for 30, 50, 100 and 200 fb~ x of integrated 
luminosity. To calculate the limits, we have used the in- 
variant mass as the discriminating variable. Since the 
distribution has a long tail, the analysis has been re- 
stricted to M 77 < 1.5 TeV, so as to have sufficient events 
for the considered luminosity. The lower limit in the M 77 
was essentially determined by the requirements on P^,. 

We have checked the stability of the limits and found 
that the 99% CL values suffers only a very small error 
(< 0.5%) from the uncertainty in the position of the Xb 
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FIG. 10: x 2 distribution for S + B and B type hypothesis for 
a given A - M q * point with 10 6 MC trials at 100 fb~ L of inte- 
grated luminosity. Here 5 + B is fitted with \ 2 distribution. 
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FIG. 11: Achievable exclusion contours in the A — M q * param- 
eter space corresponding to different integrated luminosities 
at the LHC. The regions below the curves can be ruled out 
at 99% C.L. 



peak as determined from Monte Carlo trials. To find 
the dependence on the choice of kinematical cuts, we re- 
duced the fiducial volume from \i]\ < 2.5 to \r/\ < 1.5. 
This changes the CL from 98% to 99% CL. Similarly 
the 98% CL limits obtained with P} > 200 GeV changes 
to 99% CL at P} > 250 GeV but at the cost of severe loss 
in signal efficiency. Since we have used the deviation of 
the invariant mass from the SM prediction as a discrim- 
inating variable, we expect to further improve the limit 
by combining some other uncorrelated variables [26j. 



IX. SYSTEM ATICS 

As described in the earlier sections, we have performed 
a detailed analysis including a realistic simulation of the 
various detector effects and uncertainties. Some system- 
atic uncertainties persist still and, in this section, we 
present an estimation for each of these. 

• Choice of PDF: To estimate the uncertainty due 
to the choice of the PDF, the cross sections 
were calculated with different choices of PDFs 
and the results obtained compared with those ob- 
tained for CTEQ6M (2?| ■ For comparison we used 
CTEQ5M1, CTEQ5L and MRST2001. A maxi- 
mum uncertainty of ~7% was found when CTEQ5L 
was compared to CTEQ6M. For CTEQ5M1 and 
MRST2001 these values are 2.3% and 3.5% respec- 
tively. 

• Scale Variation: To estimate this, the factorization 
scale Q (chosen to be y/1 in our analysis) was varied 
from in the range Q 2 £ [s/2, 2 s]. Also used was 
Q 2 = Pj,. In all these variations, the maximum 
uncertainty was found to be 1.6%. 

• Higher-order effects: The SM processes relevant 
to us have been studied in the literature at great 
length. Most higher order effects can be ade- 
quately parametrized in the form of a if-factor. 
For true diphoton production, these are 1.5 (Born 
process) and 1.2 (box) [3fl|. For the 7 + jet 
events, these are 1.66 when the quark fragments 
into a photon [3(3] and 1.0 when an (almost) iso- 
lated 7T° in the hadronic jet fakes a photon [301 ] . 

For the new physics contribution, the if-factor 
is not known though (indeed, the very definition 
could be ambiguous for a nonrenormalizable the- 
ory), and hence we have not used any in our anal- 
ysis. However, in the limit of a very large M q * , the 
new physics effect should be describable in terms of 
an effective operator involving quarks and photons 
and the if-factor, in this limit, is not expected to 
be too different from the SM one (T2I ]. 

If one assumes the signal if -factor to be indeed 
similar to the overall background one, then the net 
effect is a scaling of Eq.Q by a factor of K . This 
translates to a modification in the separation be- 
tween the peaks of the two histograms in FigfTUlbv 
a factor of K and is equivalent to an increase in the 
luminosity by the same factor. To be conservative, 
we choose to ignore the consequent improvements 
in the exclusion limits. 

• Energy resolution: To study the effect of the de- 
tector energy resolution on this analysis, the en- 
ergy of the photons was smeared with the stochas- 
tic term of the CMS electromagnetic calorimeter 
energy resolution [2]]]. The effect was found to be 
negligible. 
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• Dijet background: Due to limitations in computing 
resources, we did not fully simulate the background 
from jet-jet events. Although the dijet cross sec- 
tions are very large, given the low probability of 
a jet faking a photon (as described earlier in the 
text), it is obviously not very likely that too many 
such events would survive the selection criteria that 
we have imposed. A parton-level Monte Carlo cal- 
culation readily verified this. 

Even in the corresponding PYTHIA study, it was 
again observed that the kinematical and isolation 
cuts reduces this background drastically. In a sam- 
ple of 9000 jet-jet events, no event survives the 
final selection requirements. However, with the 
same survival efficiency as for 7+ jet events (i.e.,^1 
%) and with same kinematical and isolation cuts, 
we expect to have a jet-jet background of less 
than 3.7 events for an integrated luminosity of 1 
fb^ 1 . Hence we may safely assume that two photon 
events from jet-jet background will have negligible 
effect on the final confidence level calculation. 

• Luminosity error: At the LHC, for an integrated 
luminosity above 30/6 -1 , the error on the measured 
luminosity is expected to be 3% [28]. 

We have determined the effect of uncertainty in the 
theoretical cross-section on the CL. To get a con- 
servative estimate we lowered the cross section by 
1% and found that 99% CL changes to 98% CL. 

X. CONCLUSIONS 

To summarise, we have investigated the potential of 
using the diphoton final state at the LHC in probing pos- 
sible substructure of quarks. In any model of quark com- 
positeness, excited states occur naturally and these cou- 
ple to the SM counterparts through a generalised mag- 
netic transition term in an effective Lagrangian. Conse- 
quently, the presence of such states would alter the dipho- 
ton cross section, the extent of which depends on both the 
mass M q * and the compositeness scale A. The deviation 



concentrates in the large pr regime, especially for larger 
Mg* and can be substantial. For example, A = M q *=l 
TeV leads to a ~12% deviation in the cross section (when 
restricted to an appropriate part of the phase space as 
defined in Section IV). 

Using the photon reconstruction algorithm as used for 
the CMS dectector at the LHC, we perform a realistic 
estimation of the deviation caused by the excited quark 
exchange contribution to the diphoton rate. We have ac- 
counted for all major backgrounds to evaluate the limits 
in the A — M q * parameter space. The possible exclusion 
limits are very strong and depend only weakly on the 
choice of the kinematical cuts. 

While direct searches can lead to very strong limits 
from the non-observation of mass peaks, the search strat- 
egy outlined here can prove to be a complementary tool. 
In particular, as shown above, this mode is sensitive to 
excited quark masses far above the kinematical limit for 
pair-production (which mainly proceeds through gauge 
interaction). Furthermore, this method is sensitive to 
the magnetic transition coupling (q*qj) in an unambigu- 
ous manner free from all other couplings and parameters 
of this essentially complex theory. 
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